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Abstract. We present new 450 and 850 ^im SCUBA data and 1.3 mm MAMBO data of the dwarf galaxies II Zw 40, 
^ . He 2-10 and NGC 1140. Additional ISOCAM, IRAS as well as ground based data are used to construct the observed 

' mid-infrared to millimeter spectral energy distribution of these gala xies. These spe ctral energy distributions are 

, modeled in a self-consistent way, as was achieved with NGC 1569 JCalliano et alJl2003l) . synthesizing both the 

■ global stellar radiation field and the dust emission, with further constraints provided by the photoionisation of the 

gas. Our study shows that low-metallicity galaxies have very different dust properties compared to the Galaxy. 
, Our main results are: (i) a paucity of PAHs which are likely destroyed by the hard penetrating radiation field, 

a very small (~ 3 — 4nm) average size of grains, consistent with the fragmentation and erosion of dust particles 
by the numerous shocks, (in) a significant millimetre excess in the dust spectral energy distribution which can 
, be explained by the presence of ubiquitous very cold dust (T = 5 — 9 K) accounting for 40 to 80 % of the total 

PL(' dust mass, probably distributed in small clumps. We derive a range of gas-to-dust mass ratios between 300 and 

2000, larger than the Galactic values and dust-to-metals ratios of 1/30 to 1/2. The modeled dust size distributions 
are used to synthesize an extinction curve for each galaxy. The UV slopes of the extinction curves resemble that 
^ observed in some regions in the Large Magellanic Cloud. The 2175 A bumps of the modeled extinction curves are 

, weaker than that of the Galaxy, except in the case of II Zw 40 where we are unable to accurately constrain the 

2175 A bump carrier. 
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1. Introduction tinction in galaxies, in order to explore the effects of 

different environments on the dust properties. We be- 
Since dust can either reveal or conceal star formation g^^^ gg-^^^g ^^^^ metallicity environments. What 
activity, it plays an important role, often as an elusive ^een gleaned to date indicates that low metallicity 
villain, in the interpretation of the cosmic star forma- environments can, indeed harb or non-negligible q uanti- 
tion history. While characterisation of the detailed dust ties of dust llXhuan et al.lll999t IPlante fc Sauv"^E)ni 
properties is required to understand the attenuation of iGalhano et al.ll200.1^ . the effects of which can not be ig- 
the starlight, and to determine the in situ star forma- ^^^^^ Additionally, the physical pro perties of the d ust do 
tion properties, it is difficuh to achieve this due to the not resem ble those of the Galaxy JCxalliano et 
necessity of sampling a broad range of the dust spec- |Vanzi & Sauvage .27)01 . These resuhs cah to question 
tral energy distribution (SED). We have been carrying ^^e validity of the seemingly-innocuous assumption of 
out studies of detailed SED modeling, with the goal of Galactic-dust properties. AU of these resuhs have striking 
reproducing a self-consistant model for emission and ex- implications on the interpretation of extragalactic SEDs, 
and influence our view of the cosmic star formation his- 
tory. 
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NASA Goddard Space Flight Center, Greenbeh, MD 20771, This is the third paper in a series of publica- 

USA tions studying the dust properties in star-bursting low- 
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metal licity dwarf galaxies . T he first paper (hereafter pa- 
per L iMadden et al.l2005l see lMaddeiJl2000l for a preview 
of this study) was a mid-IR spectroscopic study of a sam- 
ple of dwarf galaxies obser ved by ISO CAM. The second 
paper ijGalhano et alJl200^ hereafter paper II) presents a 
detailed model of the dust SED of the dwarf NGC 1569. 
Here, we present a study similar to paper II, applied to 
the three starbursting dwarf galaxies II Zw 40, He 2-10 and 
NGC 1140. In this paper we pull together all that we have 
learned from in depth analyses of an assortment of low- 
metallicity environments. 

IIZw40 is a blue compact dwarf galaxy (BCD), at a 
distance of D ~ 10 Mpc, often considered to be the pro- 
totypical Hll galaxy, its nucleus being dominated by one 
large Hii region of ~ 0.5 kpc diameter. Its bright opti- 
cal core is located inside a large Hi envelo pe. The two 
tails seen in optical an d near- infrared (e.g. ICairos et al.l 
l200lt IVanzi et aI ."l996V arc the relics of the merger of 
two smaller galaxies (Brinks fc Klein.1988^1 . The starburst 
could be very young. The metallicity has been deter- 
mined to be Z - 1/6 Zq ((Masegosa et al.lll994l). This 
galax y is H i rich {M (H i) = 4.4 x 10^ M,7. : Ivan Zee et all 
Il998|) . while the mole cular gas emissio n is very weak 
(M^oi < 0.4 X 10^ Mq : lMeier et al1l200l|) . The radio con- 
tinuum emission is compact and mainly thermal free-free 
l|Beck et al.ll2002j) . 

He 2-10 (ES0 495-G21) is a very bright southern BCD 
galaxy, located at a distance of -D ~ 9 Mpc. It was 
the fi rst Wolf-Rayet galaxy to be identified ijAllen et al.l 
Il976l) . Its peculiarities among BCD g alaxies is that its 
meta llicity appears to be almost solar ijKobulnickv et alJ 
l1999^ an d it has a relative ly large amount of molec- 
ular gas l)Meier et alJ l200l|) . The starburst is concen- 
trated in two regions separated by a dust lane. It con- 
tains several Su per Star Clusters (SS Cs) probing a burst 
of a few Myr l).Tohnson et al.l lioool) . This recent burst 
could have been triggered by the inte raction with a 
molecular cloud llKobuliiickv et al.lll995(l . A radio map 
i Kobulnickv fc JohnsonI Il999(l reveals several compact 
thermal continuum sources, indicating ultra-dense Hii re- 
gions. Ha images sho w a large scale bipolar outflow of gas 

NGC 1140 is an amorphous starburst galaxy, larger 
than IIZw40 and He 2-10 with a larger reser voir of gas. 
Its size is roughly 1.7' x 0.9' l)Buat et alJl2002l) . at the dis- 
tance of D ~ 20 — 25 Mpc. The central Hil complex of 
NGC 1140 has an Ha lumi nosity which is 40 tim es higher 
than that of SODoradus I Hunter et al.lll994a^ and con- 
tains several SSCs w hich have been stud i ed in detail by 
iHunter et aP l|l994hl) and Ide Griis et alJ l)2004fl . Its star 
formation rate is about 1 M(?,yi~^. NGC 1140 may be in 
the final stages of a merger ( Hunter et alJl994a|) and has a 
small dwarf companion. Among the four dwarf galaxies of 
our sample (NGC 1569, IIZw40, He2-10 and NGC 1140), 
it is the only one which shows pronounced PAH emission 
at mid-IR wavelengths (paper I) . 

The modelling of the global SED of a dwarf galaxy, 
using actual dust models, has been presented bv iMaddeni 



(|2000l), ' Galhano et all ll2002^. iLisenfeld et all 
ICalliano et alJ tOo£^ iTakeuchi et alJ J2003irand this 
study (see also iGallianol l2005l) . iTakeuchi et all l(200.1^ 

have used a d ust evolution inqdel t o fit the IR SED of 
SBS 335-052. iT.i.senfeld et tOO'j) and lGa,1lia,no et a,lJ 
l|2003|) (pa per II) applied the iDesert et all l)l990|) model to 
NGC 1569. iGalhano etal] l)2003(l . having sufficent obser- 
vational constraints, were able to constrain the dust size 
distribution of different dust components for NGC 1569. 
Our conclusion in paper II was that the global IR SED 
in NGC 1569 is composed of dust components with very 
different size distributions from our Galaxy and con- 
trary to our Galaxy, is dominated by small dust particles 
(< 10 nm). Additionally, a submillimetre excess was found 
that could be attributed to a very cold dust component. 

Studies on the ISM of dwarf gal axies, on smaller spa- 
tial scales have also been conducted. iBot et al.l ll2004l) ap- 
plied the iDesert et alJ l)l990l) model to the diffuse emis- 
sion of the Small Magellanic Cloud (SMC). While they 
did not vary the size distribution, they took into account 
the variations of the radiation field. They had di f ficulty 
fitti ng the 60 Mm flux . Finally, IPlante fc Sauvagd l)2002() 
and lVanzi fc Sauvag^ |2004) have modelled the dust SED 
around the SSCs in SBS 0335-0 52 and NGC 5253 respec- 
tively. They used Dusty (Jvezic fc Elitziirl Il997l) which 
solves the radiative transfer equations in a spherical en- 
vironment, without considering the process of stochastic 
heating. In these cases, the size distributions of the local 
environments were found to be dominated by large grains. 
These differerent results obtained for global galactic scales 
versus more local regions around the SSCs is likely due to 
the destruction of small grains around SSCs, as is also o b- 
served in the vicinity of AGNs (|Maiolino et a l."2001a^lJ). 

The paper is organised as follows. Section|21presents an 
overview of our new observations and the data we obtained 
from the literature. Section 0] presents the modeled SEDs 
of II Zw40, He 2-10 and NGC 1140 and the consequences of 
the results on the dust properties. We end with a summary 
and the conclusions in Sect.|Sl which ties together Papers 

1, II and III. 

2. The observations 

Observed SEDs were constructed as completely as pos- 
sible, using our ISOCAM data (paper I), our new 450 
and 850 /im JCMT (SCUBA) observations and our IRAM 
(MAMBO) observations presented here, as well as incor- 
porating data from the literature for various telescopes 
and wavelengths. The differences in beam sizes are not of 
great concern for this paper, since we are modeling the 
global SEDs here. 



2.1. SCUBA images 

We obtained 450 and 850 ^m data of IIZw40, H e 2- 
10 and NGC 1140, with SCUBA llHofland et al.lll999|) . a 
bolometer array on the James Clerk Maxwell Telescope 
(JCMT), during two observing runs during February 2000 
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II Zw 40 (SCUBA 850 juni) II Zw 40 (SCUBA 450 juni) 




02^ 54™ 36^ 34^ 33^^ 02^ 54™ 36^ 34" 32^ 



Fig. 1. SCUBA maps of the throe galaxies. The first row of images is IIZw40, the second is He 2-10. and the third 
is NGC 1140. The wavelength of the left column images is 850 fxm and the wavelength of the right column images is 
450 yum. The field of view is the same for the 450 /xm and 850 /um images of each galaxy and the color table has the 
same dynamic range. The contours are 99, 90, 70, 50, 40 and 30 % of the peak flux values. 
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and December 2000. Observations were carried out in 
the jiggle-mapping mode using a 64-point jiggle pattern 
with a chop throw of 150". The precise data reduction 
method used is described in detail in paper II includ- 
ing details of the error analyses. The full width at half 
maximum (FWHM) of the beam is 8.5" at 450 /im and 
15.2" at 850 /ini. The measured tcso at 225 GHz, obtained 
from the Caltech Submillimeter Observatory (CSO) ra- 
diometer, ranged from 0.04 to 0.1 during our observa- 
tions and our observed calibration sources were Uranus, 
Mars or CRL 618. The final images are shown in Fig. ^ 
and the fluxes are presented in Table ^ The signal-to- 
noise ratios of these images are: S'/A^(850 /im) ~ 10 and 
S/N{i50iim) ~ 7, for NGC1140; S/N{850fim) ~ 15 and 
S/NlA50fim) ~ 11, for IIZw40; S/N{850fim) ~ 23 and 
S/N{A50iim) - 14, for He 2-10. 

The contamination of the 850 fim fluxes, due to the 
CO (3-2) line, can be quantified from observations re- 
ported in the literature. Based on CSO observations, 
[Meie r ct al. ( 200^ give an upper limit for II Zw 40 of 
Jco(3-2) S 0-9 Kkms~^ in 22", which we convert to 
^C^O{3^) ^ 3.4 mJy, for the total galaxy by using the de- 
convolved 850 /im SCUBA map and assuming that the 
CO emissi on has the same spatial distribution as the 
cold dust. iMeier et"al ] ll200lh also give the /co(3-2) for 
several pointings in He 2-10 which have some overlap 
with our SCUBA data. We took the central measurement 
/co(3-2) = 16-6 Kkms~^ and scaled it in the same man- 
ner as for IIZw40. The estimated CO(3-2) flux in the 
850 /xm band for the total galaxy is, thus, -F'co(3-2) ^ 
30 mJy (TableHJ. In the case of NGC 1140, we did not find 
any CO(3-2) measurement, in the literature. However, this 
contribution should be very low considering the contribu- 
tion of the C0(2-l) line in the MAMBO band fSect. l^ . 
thus we believe that this contribution is not significant. 

The radio continuum contribution to our total sub- 
millimetre fluxes can be estimated by extrapolating from 
the numerous radio fluxes given in the literature. The ra- 
dio continuum at submillimetre wavelengths is normally 
domi nated by free-free emission in galaxies (e.g. ICondonI 
Il992tl . The free-free emis sion is described by F , j oc 
In the case of IIZw40, ISramek fc Weedmai] lll986|) re- 
port a global thermal flux of Si, = 13.9 — 15.2 mJy at 
4.8 GHz. We deduce a flux S"^ = 9.1 - 9.9 mJy at 850 //m 
(< 10%) a nd S„ ^ 8.5 - 9.3 mJy at 45 um (< 4%). 
For He 2-10, iKobiilnickv fc .TohnsonI l(l999^ report a total 
flux of Fi, = 21.1 ± 1.2 mJy at 14.9 GHz, which we con- 
vert to S't. = 14.5 - 15.4 mJy at 850 /xm (< 13 %) and S^, = 
13.6- 14.4 mJy (< 4%) at 450 /im. We could not find a de- 
composition between thermal and non-thermal flux, thus 
this estimate serves at the upper limit of the radio contri- 
bution in our bands. We remain conservative here in order 
to be sure that if we do observe a submillimetre excess 
arising from dust in these galaxies, as found previously in 
NGC 1569 (paper II), the results will not be effected by 
our assumptions of contamination from non-dust sources. 
Finally. iKlein et al.. ((1983.) report a = 11 ± 3mJy flux 



at 10.7 GHz, in NGC 1140. Thus, we are able to give a 
range of values of 5^ = 5.6 — 9.9 mJy at 850 //m (< 12 %) 
and = 4.9 - 8.6 mJy at 450 /xm (< 4 %). These contri- 
butions are summarised in the last part of Tabled Thus, 
subtracting the radio continuum and CO contributions to 
the observed submillimetre fluxes, we are confident that 
we can quantify the flux density arising from dust only. 





IIZw40 


He 2-10 


NGC 1140 




(mJy) 


(mJy) 


(mJy) 


77'tot 

^450 fim 


248 ± 81 


342 ± 65 


272 ± 55 


TTltot 

-f^850 fim 


98± 14 


130 ± 12 


69 ±28 


CO(3-2)850 /im 


^ 3.4 


~ 30 




Radio: 








- 450 fim 


8.5-9.3 


13.6 - 14.4 


4.9 - 8.6 


- 850 fim 


9.1 - 9.9 


14.5 - 15.4 


5.6-9.9 



Table 1. Flux and contributions to the submillimeter 
bands of II Zw 40, He 2-10 and NGC 1140. The total fluxes 
are not corrected for radio continuum and CO (3-2). 



As noted for NGC 1569 (paper II), the morphology of 
the submillimetre emission (Fig. ^| is not concentrated 
toward the outer regions. It follows the mid-IR emission 
relatively well, the peak correlating with the star forming 
regions. 

2.2. MAMBO observations 

We obtained 1.3 mm data (230 GHz, bandwidth~80 GHz) 
of IIZw40. He 2-1 and NGC 1140, with MAMBO 
ijKrevsa et al.l Il99^ . a bolometer array on the IRAM 
30 meter radio-telescope, during two observing runs in 
December 2001 and January 2002. Observations were car- 
ried out in On-Off mode using the 117-bolometer array. 

The On-Off observations were conducted in standard 
chop-nod mode, with individual scans divided in 6 sub- 
scans each of which yielded 20 seconds of on-|-off source 
exposure. The secondary mirror was chopped by 55" in 
azimuth for IIZw40 and He 2-10 and 70" for NGC 1140. 
The pointing and the focus of the telescope were checked 
every hour. The data were analysed with the NIC soft- 
ware llBroguiere et al.ll2002|) from the GILDAS package. 
The data reduction included the following steps. 

1) Atmospheric opacity: The skydip measurements were 
carried out every two hours. The T230, optical depth of the 
atmosphere at 230 GHz, varied between 0.14 and 0.35 dur- 
ing the observations and was below 0.2 during most of the 
scans. Moreover, it was relatively constant between suc- 
cessive sky dips. Each scan was corrected for atmospheric 
extinction using linear interpolation of the opacity mea- 
surements obtained from the skydip preceding the scan 
and that following the scan. 

2) Gain and bad bolometers: Some bolometers are known 
to be noisy. We flagged these currently-identified bad 
bolometers during the remainning procedures. 

3) Spikes: A despiking function is applied to each bolome- 
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ter, with a threshold of 5 x cr. The spikes found are removed 
and replaced by the interpolation of the adjoining data. 

4) Final signal: The mean signal of each bolometer and its 
variance is computed from the corresponding subscans. A 
baseline is removed from the raw signal by a least-square 
fit of a straight line to the weighted sequence of ON and 
OFF means. The high-frequency noise is removed from all 
of the channels. Finally, the flux of the reference channel 
is determined by removing the flux density of the adjacent 
four channels. If Tib is the flux in counts of the scan i in the 
bolometer h, we calculate the mean, Tb = {Tib) scans, of ev- 
ery scan flux, Tb, weighted by 1/crf , where ai — <j{J-'ib)hoio 
is the deviation of each scan. We subtract the foreground 
emission which is determined from the mean flux of the 
bolometers which are not on the source. Thus, the flnal sig- 
nal of each bolometer is Tl;°" = Tb — {Tb)i,__^^^^j.^^y These 
j^corr plotted in Fig.|21for the two galaxies which are 
detected. 

5) Calibration: Mars, a Ori, NGC 7538 and HL r were ob- 
served daily as primary and secondary calibrators. The 
calibration data were reduced in the same way as other 
galaxies. We compute a flux conversion factor, $, which 
provides the relationship between counts and astrophysi- 
cal fluxes (in Jy), averaged over the calibrators. The flux 
of our sources is converted to astrophysical flux values 
(in Jy/beam), = $ x T^°" , by multiplying the signal 
by the conversion factor. We use the submillimeter spa- 
tial distribution to scale this flux to the total flux of each 
galaxies. First, we deconvolve the SCUBA 850 /Ltm maps 
with a n iulti-res olution Lucy algorithm from the MR/1 
package l)Starck et al. 19981 . using a 3o- detection thresh- 
old. Second, we convolve this image to the IRAM 30-meter 
beam and we extract the ratio of the total flux to the flux 
in one beam, TZ = F*°*(850^m)/F|-,oam(iRAM)(850^m). Wc 
multiply the IRAM flux in Jy/beam by this ratio to ob- 
tain the total flux at 1.3mm, F*°*(1.3mm) = Fb^^^ x TZ, 
where 6ref is the reference bolometer. The final fluxes are 
given in Table |21 the values of TZ is given in the last part 
of this table. 

6) Error estimation: The error on the total flux is given 
by AF*°*(1.3mm) = AF^,„J7^ + Fb^^^ATZ. The error on 
the aperture, ATZ, is estimated by shifting the SCUBA 
image convolved by the IRAM beam by ±3" in RA 
and DEC. This 3" error is the typical IRAM point- 
ing error llWildlll999 ') and is consistent with deviations 
observed during the pointing verification measurements 
during the observations. The error on the flux in one 
beam, AFb^^^ , is the quadratic sum of two components: 
(AFb^J^ =V(F,,_|,,„,e}) + (A$/$)^ The flrst compo- 
nent, cr(-Fb-{sourcc})j is the RMS of the final signal. The 
second component, A$/$, is the relative deviation of the 
flux conversion factor over the various calibrators. The 
various contributions to the error are given in percentage 
in Tabled 

7) Non-dust contamination: The expected contamination 
from the CO (2-1) line can be quantified from measure- 
ments in the literature. In the case of IIZw40, we esti- 
mate 0.1 mJy due to the C0(2-l) fine in the IRAM 1.3 mm 




20 40 



Channel 
He 2-10 




Fig. 2. On-Off observations for the two detected galaxies, 
IIZw40 and He 2-10. These plots show the signal in each 
bolometer after foreground subtraction. The peak corre- 
sponds to the reference channel (n° 53) which is centered 
on the source. The other channels are assumed to see the 
sky. 



broadband from ISaee et al ] lll992l) . For He 2- 10, we esti- 
mate a flux of 4mJy in a 12.5" aperture from lBaas et alJ 
l|l994|) , and scale it to ^ 5 mJy for the total galaxy. The 
upper limit given bv lHunter fc Saed il993), for NGC 1140, 
gives a C0(2-l) flux in the IRAM band of < 0.02 mJy. The 
radio continuum component can be extrapolated from ra- 
dio fiuxes that we used in Sect. 12. ll In the case of II Zw 40, 
we deduce a fiux 5^ ~ 9.5-10.3mJy at 230 GHz (< 24%); 
for He 2-10, < 15.1 - 17.0 mJy at 230 GHz (< 28%) and 
~ 5.9 - 10.3 mJy at 230 GHz for NGC 1140 (~ 20%). 
These contributions are summarised in Table |2] 

2.3. Infrared data 

The details and the data treatment of the ISOCAM CVF 
mid-IR observations are presented in paper I. We found 
that the slope of the mid-IR spectrum is a critical factor 
in constraining the dust model, particularly in the 5 to 
16 /zm wavelength range (paper II). Thus, using only the 
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IIZw40 


He 2-10 


NGC1140 




(m ^v) 


(m ^v^ 


(injy) 


F*°*(1.3mm) 


43 ± 13 


60 ± 14 


< 48 


Radio cont. 


9.5 - 10.3 


15.1 - 17.0 


5.9 - 10.3 


00(2-1) line 


0.1 


4 


< 0.02 


RMS error 


17% 


2% 


100% 


A$/$ 


16% 


12% 


16% 


A7^ 


22% 


25% 


30% 


l/7^ 


41% 


48% 


22% 



Table 2. Contributions to the 1.3 millimeter fluxes 
of IIZw40, He2-10 and NGC 1140. The total fluxes, 
i^^°*(1.3 mm), are not corrected for radio continuum and 
C0(2-l). 



IRAS 12 fj,m band does not provide sufficient constraints 
on the mid-IR SED which predominantly traces the hot 
dust component. From a fit of the spectrum in II Zw 40, we 
find that the internal extinction is Ay ~ 15 (paper I). We 
use the deredenned spectrum in the modeling. We char- 
acterised the mid-IR dust continuum by choosing a few 
wavelength regions of the CVF which do not contain aro- 
matic bands or ionic lines (Table 0. 

IRAS fluxes for II Zw40 ar id NGC 1140 ar e give n 
bv iHunter et all lll989t) an d by iMelisse fc Israeli l|l994j) . 
iTh ronson fc Telescol l|l986|) also give the IRAS fluxes for 
IIZw40. The values for the four IRAS broadbands given 
by these authors are consistent. Moreover, the 12 /xm 
IRAS flux is consistent with the ISOCAM CVF spec- 
trum integrated into this band (pa per I). For ou r mod- 
eling purposes, we use the values of iHunter et all l|l989() 
since error bars are also provided. These fluxes have been 
color-corrected. We u se the IRAS fluxes for He 2-10 from 
ISauvaee et alJ l|l997|) . 

In the case of II Zw 40, an ad ditional 20 /im flux of 
1.0±0.2Jy is given bv iRoche et all (1991). For He 2-10, we 
consid er the flux at 11.65 /^m reported bv fSauvaee et al.1 
who also give a N band flux for He 2-10. However 
we use the N band flux from lVacca et alJ (^002) instead, 
since there is a large discrepancy between this flux and 
that cited bvlSauvage et aD l|l997|) (see the discussion by 
IVacca et al.l2002() . We include an unpublishe d M band ob- 
servation of He 2-10 ((Sauvage fc Vanzill2005l) in our SED. 



2.4. Optical data 

To constrain the input stellar radiation field we use op- 
tical data from the literature. There are several papers 
reporting optical observations of IIZw40, He 2-10 and 
NGC 1140. We choose magnitudes explicitely given for the 
total galaxy and where the authors specified the Galactic 
foreground extinction corrections. Preference was given 
to measurements quantifying the uncertainties. In some 
cases, we were obliged to scale data which were given for 
an aperture smaller than the total size of the galaxy. 

In the ca se of II Zw 40 , we u sed B, V, R and I mag- 
nitudes from ICairos et al.1 l)200lh and U magnitude from 



the RGB catalog llde Vaucouleurs e~ [Jl99lh . The B mag- 
nitude, given bv ICairos et alJ ll2001 ) is als o consistent 
with t hat p rovided bvlHeisler fc Vadcr (' l994^ ■ lDeeg et all 
ll1997^ and EI Vaucouleurs et al.. ((1991.) . The V magni- 
tude of ICairos et alJ ll200lh is consistent with that of 
Ide Vaucouleurs et al.l lll99ll) an d the R and I magnitudes 
are consistent with iDeeg et all lll997l). Th e J, H and K 
fluxes of II Zw 40 reported by Vanzi et alJ ()l996(l for an 
aperture of 27", were scaled to obtain values for the entire 
galaxy. We computed the scaling factors from the original 
images l)Vanzi et al.lll996|) . 

He 2-10 has been observed bv IJohanssonI ()l987|) in U, 
B, V, R, I, J, H and K bands in several apertures. U, B 
and V observations are given for the total galaxy (aper- 
ture of 61"), while the other bands have been observed 
in, at most, a 31" aperture. The difference is relatively 
small (15 %). We scaled the 31" flux measurements to fit 
the U, B, V bands in the total galaxy aperture. The error 
bars were not provided by the author. Consequently, we 
took the deviation of the two sets of observations reported 
by the author and added the error induced by the differ- 
ences of apertures. Even with these considerations, the 
uncertainties we obtained are relatively small and may be 
underestimated. 

F or NGC 1140, we used U, B , V, R and I magnitudes 
from iGallagher fc Hunted lll987t) and J , H and K magni- 
tudes from lHunter fc GaflagheJ i)l985|) . iBuat et all l|2002ll 
also observed NGC 1140 around 0.1 /zm, in a 30" x 30" 
aperture which we believe encompasses the majority of 
the emission of the galaxy at these wavel engths, which 
is conc entrated within a diameter < 3 kpc l)Hunter et alJ 
Il994bl) . 

These data, summarised in Table are used to con- 
strain the stellar SED. 





IIZw40 


He 2-10 


NGO 1140 


Band (/im) 


Flux (mjy) 


Flux (mJy) 


Flux (mJy) 


0.36 (U) 


15.3 ± 2.3 


63.2 ±2.1 


15.1 ± 2.3 


0.44 (B) 


23.8 ± 1.6 


90.8 ±3.0 


21.2 ±3.2 


0.55 (V) 


32.9 ±6.3 


100.7 ±5.0 


23.6 ±3.5 


0.64 (R) 


36± 11 


122.1 ±4.0 


26.8 ±4.0 


0.79 (I) 


29.3 ±9.7 


126.5 ±4.1 


35.2 ±5.3 


1.26 (J) 


20.6 ±4.1 


147.9 ± 4.8 


36.4 ±5.5 


1.60 (H) 


21.9 ±3.3 


175.4 ± 5.7 


40.1 ±6.0 


2.22 (K) 


18.5 ± 2.8 


130.2 ± 7.5 


28.6 ±4.3 



Table 3. Optical data from the literature (references in 
the text). These fiuxes are global values, corrected for the 
foreground Galactic extinction but not for internal extinc- 
tion. 



3. Self-consistent modeling of the global SED 

The modeling of the global SEDs of II Zw40, He 2-10 and 
NGC 1140 is done using the exact same procedure as de- 
scribed in paper II, for the case of NGC 1569. We refer to 
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this paper for the detailed description of the procedure we 
use here. The main steps are the following. 

— UV-to-optical data from the literature (Table EI), cor- 
rected for the Galactic extinction are used to constrain 
the stellar radiation field. This radiation field is mod - 
eled with PEGASE (Eioc fc Rocca-Volmerang3ll997|) . 
a stellar evolutionary synthesis model. We consider two 
instantaneous bursts of star formation representing the 
young stellar population created by the recent star- 
burst and the old underlying population. We are only 
interested in the shape of the ISRF here, and the age is 
a free parameter. Hence, the initial metallicity is not 
of great concern, due to the ag e-metallicity degener- 
acy (e.g. iLe Borgne et al horn . However, for a given 
metallicity, several age combinations fit the observa- 
tions. 

~ The stellar population age degeneracy is removed 
by constraini ng a photoionisation model (CLOUDY; 
lFerlandlll99l with the mid-IR ionic line ratios that 
we measure (paper I). 

— This radiation field is used to hea t the d ust which 
is modeled with the IPesert et all l|l99(l hereafter 
DBP90) model. We fit the observed IR SED by vary- 
ing the dust size distribution. The solution allows us 
to synthesize an extinction curve self-consistent with 
the emission properties. 

— We iterate this process, correcting the UV-to-optical 
data for internal extinction with the synthesized ex- 
tinction curve, until we reach an agreement between 
the extinction and emission. Only the spectral depen- 
dency of the extinction curve is taken into account. 
The global internal optical depth is deduced from the 
energy balance between stars and dust: 



UV-opt ■ 



(1) 



where F^, is the intrinsic flux emitted by the stars, 
^uv-opt, the escaping stellar flux, and i^iR_mm, the flux 
reemitted by the dust (these fluxes are integrated over 
frequency) . 

3.1. The observed SEDs 

The stellar SEDs are constrained by the fluxes in Table 13 
and by mid-IR ionic line ratios (paper I). The [Neii] 
(12.81 /Ltm) line is not detected in IIZw40, thus, in this 
case, we constrain the photoionisation model with only 
two ratios: [S iv]/[Ne ill] and [Ariii]/[Neiii]. In the case 
of He 2-10, w e do not have a n ISOCAM CVF spec- 
trum. However lBeck et al.Ml997j) report NASA IRTF mid- 
IR spectroscopy of this galaxy and they measure ionic 
line intensities. We deduce from this study, the ratios 
[Neil] /[Arm] = 12.5 ± 5.2 and [Ncii]/[Siv] = 50 ± 40. 

The dust SEDs are constrained by the data points 
shown in Table 01 All the non-dust contributions have al- 
ready been subtracted from these fluxes listed in Table 0] 



3.2. The dust model 

As in paper II, we define an effective radius of the galax- 
ies, ReS, which is the effective distance between stars and 
dust in a thin shell model. We estimate the equivalent ra- 
dius, i?equiv, by fitting an ellipse to the 850 /im images and 

ReS = 3/4 X i?equiv 

For II Zw 40, we assume a distance of -D = 10 ± 1 Mpc 
and for He 2-10, we assume D = 9 ± 1 Mpc, consis- 
tent with all the rece i it values adopted in the litera- 
ture llRoche et al.lll99lUContilll99ll: FVacca fc Contlll992l: 
iTelesco et al.lll993|) . In the case of NGC 1140, we adopt 
D = 23 ± 2 Mpc to be consistent with the recent values 
of 20 to 2 5 Mpc ("Kinn ev et alJll 993: Hi mter et al.lll994jj 
ISaikia et aL.1994.Heckman et al..l998;.Buat et alJl2f 
The relevant model parameters for the sources are given 
in Table 



Name 



R^B D M(Hi) Z/Zc^ 

(kpc) (Mpc) (Mq) 



II Zw 40 
He 2-10 
NGC 1140 



0.5 ±0.2 10 ±1 
0.8 ±0.3 9±1 
2.6 ± 0.9 23 ±2 



4.4 X 10** " 
3.1 X 10* " 

7.5 X 10^ " 



1/1 
1/3^ 



Table 5. Physical parameters of IIZw40, He 2-10 and 
NGC 1140. Rcff is deduced from our mid-IR and sub- 
millimetre images, the other parameters are found in 
the literature. D is the distanc e to the galaxy an d 



Z, the metal abunda n ces, (a ) Ivan Zee et al.[ 
(6) e.g. IWalshfc Rovl |l99^ iMasegosa et al.l 



re^lHunter et a,1..(.1994 



(c) Igauvage et al.' 
)|Hi 



(d) [Kobulnickv et al.' 
; rnlGuseva et al. (,200(>'). 



1998) 



1994) 



1S22D 



The DBP90 dust model was originally designed to 
explain the dust properties of the Milky Way and in- 
cludes three dust components: the Polycyclic Aromatic 
Hydrocarbons (PAHs) which are responsible for the mid- 
IR emission features at 3.3, 6.2, 7.7, 8.6 and 11.3 /im, the 
Very Small Grains (VSGs), which are assumed to be the 
carriers of the 2175 A extinction bump and the Big Grains 
(BGs), which have an emission peak in the far-IR. VSGs 
are carbonaceous grains and BGs are silicate grains. Each 
of the 3 dust components is described by their minimum 
and maximum grain size (a_ and a+), and a, the index of 
the power law size distribution which is n{a) oc a^" (a is 
the grain radius and n(a) is the number density of grains 
between a and a + da) . The dust mass abundance of each 
component is Y = Mx/Mu, where Mx is the mass of the 
dust component, X, and Mh is the hydrogen mass in the 
galaxy. 

4. Results and discussion 

4.1. The model parameters 

We attempt to fit the observed SED with the DBP90 dust 
model, varying the dust parameters for each component. 
We retain the original denomination of BGs, even if the 
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IIZw40 




He 2-10 




NGC1140 


Instrument 


A {fJ,m) 


Flux (mJy) 


A {nm) 


Flux (mJy) 


A if^m) 


Flux (mJy) 


ISOCAM 


8.8 


89 ±55 






6.1 


41 ±26 




9.4 


132 ±47 






6.2 


89 ±22 




10.1 


177 ± 64 






6.5 


46 ±25 




10.8 


238 ± 64 






6.9 


55 ±21 




11.8 


315 ± 65 






7.4 


82 ±21 




12.1 


339 ± 67 






7.6 


119 ±20 




13.2 


417 ± 73 






7.8 


113 ±20 




13.9 


458 ± 75 






7.9 


120 ±21 




14.6 


489 ± 69 






8.1 


65 ±20 




15.0 


504 ± 78 






10.9 


82 ± 15 




16.0 


532 ± 87 






11.1 
11.3 
11.5 
11.7 
13.4 
13.9 
14.4 
16.1 


149 ± 17 
192 ± 17 
118 ± 18 

90 ±16 

91 ±23 
62 ±20 
89 ±23 
78 ±35 


Others 


20 (Q) 


1000 ± 200 


4.8 (M) 
10.1 (N) 
11.65 


60 ±9 
720 ± 95 
850 ± 100 






IRAS 






12 


1100 ± 160 








25 


2170 ± 300 


25 


6550 ± 980 


25 


390 ± 70 




60 


7280 ± 900 


60 


23800 ± 3600 


60 


3940 ± 600 




100 


6360 ± 1100 


100 


31800 ± 4800 


100 


5000 ± 800 


SCUBA 


443 


241 ± 81 


443 


328 ± 67 


443 


270 ± 110 




863 


87± 14 


863 


85 ±36 


863 


69 ±28 


MAMBO 


1200 


33± 13 


1200 


40 ± 18 


1200 


;$ 42 



Table 4. Observed dust SED. The data are the flux values used to constrain the dust model. They are corrected for 
all non-dust effects (molecular lines, radio continuum) and extinction. The references are given in the text in Sect.|21 
For NGC 1140, the ISOCAM data which are reported are those used to constrain both the continuum and the PAH 
bands. The VSGs have been constrained separately (see explanation in text). 



model results indicate that these grains are small, as in 
the cases reported here. The values of the parameters of 
the best fits to the SEDs (Fig. ^ are given in Table |B1 

The bulk of the emission originates essentially in the 
VSGs and BGs (Fig.OJ. The free parameters are a+, a and 
Y, for the VSGs, and a_, a and Y, for the BGs. We fix 
to the Galactic value, since its determination from the 
fit is very uncertain, the large sizes not being significant 
contributers to the emission. For He 2-10 and NGC 1140, 
we also fix aX^*^ to the Galactic value, since the best 
is obtained for very low unphysical values (a few A). On 
the other hand, for nZw40, the value of a_ can be de- 
termined. Among our sources. He 2-10 is the most poorly 
constrained, we do not have a CVF spectrum. However 
the number of free parameters remains smaller than the 
number of observations (TableEJ- Since the PAH emission 
is detected only in NGC 1140, Ypah remains a free param- 
eter. For IIZw40 and He 2-10, we are able to put only an 
upper limit on the PAH mass. 

Our results show that the three component DBP90 
dust model is not adequate to explain the full dust SED 
out to millimetre wavelengths for any of the galaxies, as we 
first discovered in NGC 1569 (see paper II for this point). 



We are compelled to add a fourth component for each 
of the galaxies. The Very Cold Grains (VCGs) are mod- 
eled using a modified black body, with emissivity indices 
of /? = 1 (consistent with carbonaceous grain properties) 
and (3 — 2 (consistent with silicate grain properties) . This 
very cold dust component is constrained, only by the three 
submillimetre observations. Thus it is very difficult to con- 
strain the temperature, the mass and the (3 of this compo- 
nent. We choose to vary only the temperature T and the 
mass of this component, considering the different values 
of /?. However, for II Zw 40 and He 2-10, the fit with /3 = 2 
gives a dust mass higher than the metal mass of the galaxy, 
thus we exclude this value and fit only a /? = 1 black body. 
In the case of NGC 1140, we only have an upper limit on 
the flux at 1.3 mm, so we do not have sufficient data to 
constrain the value of (3. 

Table El contains the values of several statistical quan- 
tities: n is the number of observations, m, the number of 
free parameters used for the standard model, TOvcgi the 
number of free parameters after the introduction of the 
VCGs, = X^/ {n — m — 1) is the reduced chi square, for 
the standard model, and x^vcg' the reduced chi square 
with VCGs. 
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Milky Way 


NGC 1569 


II Zw40 


He 2-10 


NGC 1140 


PAHs 


Y 


4.3 X 10"" 


< 1 X 10-" 


< 5 X 10" ' 


< 3 X 10-" 


8.3 X 10-' 




a- 


0.4 nm 


0.4 nm 


0.4 nm 


0.4 nm 


0.4 nm 




a+ 


1.2 nm 


1.2 nm 


1.2 nm 


1.2 nm 


1.2 nm 




a 


3.0 


3.0 


3.0 


3.0 


3.0 


VSGs 


Y 


4.7 X 10"" 


1.8 X IQ-'' 


2.3 X 10"'' 


4.9 X 10-" 


2.2 X 10-" 




o_ 


1.2 nm 


1.2 nm 


2.0 nm 


1.2 nm 


1.2 nm 




a+ 


15 nm 


7.8 nm 


2.9 nm 


2.6 nm 


3.0 nm 




a 


2.6 


4.0 


1.0 


5.5 


7.9 


BGs 


Y 


6.4 X 10"^ 


4.4 X IQ-* 


2.6 X 10"" 


8.9 X 10-* 


6.1 X IQ-* 




a- 


15 nm 


2.2 nm 


3.2 nm 


2.5 nm 


3.8 nm 




a+ 


llOnm 


110 nm 


llOnm 


110 nm 


110 nm 




a 


2.9 


6.3 


25 


8.9 


30 


VCGs 


Y 


unknown 


1.3 X 10"'' 


1.3 X 10"^ 


1.5 X 10-"* 


~ 1.1 X 10-^ 




T 


unknown 


5K 


8K 


6K 


~ 5K 




P 


unknown 


1 


1 


1 


~ 1 


Statistics 


n 


_ 


18 


18 


10 


24 




m 


_ 


7 


7 


6 


7 




mvcG 


- 


9 


9 


8 


10 








1.7 


4.5 


9.8 


1.5 




X VCG 




0.27 


0.36 


2.8 


0.89 




AAIC 




46 


67 


25 


52 




AAICC 




61 


82 


149 


58 



Table 6. Dust parameter values and statistical quantities, for the best model fit. For comparison, we give the corre- 
sponding values for the Milky Way (DBP90) and for NGC 1569 (paper II). 



The introduction of the VCGs, having important as- 
trophysical consequences, should be examined with scruti- 
nity. Hence, we have calculated the Akaike's Information 
Criterion (AIC), in order to examine the validity of in- 
troducing this new component. The AIC is a statistical 
quantity that can be regarded as the extension of the log 
likelihood, including the effect of the instability induced 
by the number of paramete rs. An elementa ry derivation 
of the AIC can be found in iTakeuchil tOOd). and a brief 
summary of its use in Appendix of iTakeuchi et a h* ('2000') . 
The squared sum of the deviations of the model from the 
data is: 

n 

5 = 5][log(L.,)-log(/(A.,0))]^ (2) 

being the observed monochromatic luminosity at the 
wavelength Ai, /(Ai, 9), the value of the model at the same 
wavelength, and 9, the parameter vector. The AIC is de- 
fined by: 

AIC = nx ln(27r + 1) + n x In ^-^ + 2(m + 1). (3) 

The second order corrected AICC is: 

AICC = AIC+'-^^^^±^}^^. (4) 

71 — m — 2 

A variation of AIC and AICC greater than unity, after in- 
troducing the VCGs, is required to justify this new compo- 
nent. Table|Slcontains the values of AAIC — v4/Catandard — 
AICycg and AAICC which is defined identically. We see 
that these variations are much larger than unity, thus the 
VCG component is statistically reasonable. 



The range of reliability of the parameters is provided 
in Table|7| The uncertainties of the parameters in Table[3 
are derived from both the uncertainty of the geometry 
(the value of -Rcff) and the errors of the observations 
(Section[2J). A complete description of the way we compute 
these errors is given in paper II (Sect. 4.1.2). 

4.2. The size distributions 

The size distributions are shown in Fig. 0] Compared 
to the Galaxy, there is an obvious dearth of PAHs and 
the mass spectrum is dominated by grains of small sizes 
(radius a ~ 3 — 4 nm) as shown in Fig. 01 The power- 
law size distributions for all of the galaxies are quali- 
tatively similar. The preference for grains of small sizes 
may be due to the high supernovae rate in these galax- 
ies. Supernovae produce shocks that fragme nt and erode 
the l arge grains and produce smaller grains l) Jones et alJ 
as we de monstrate in pape r II where we use the 
shock model of I Jones et all l|l996|) to explain the shape 
of the size distribution of NGC 1569. The effect is similar 
here for IIZw40, He 2-10 and NGC 1140. 

4.3. The modeled SEDs 
4.3.1. The ISRFs 

Fig-Elshows the synthesized radiation fields for our galaxy 
sample and we compare these with that of the Galaxy in 
Fig. |B| In the dwarf galaxies with evidence for current 
starburst activity, the effective, global radiation field as 
seen by the dust, is more intense and harder than that 
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NGC 1569 


IIZw40 


He 2-10 


NGC 1140 




PAHs 


Y 


~ - IQ-'' 


~ - 5 X 10" ' 


~ - 3 X 10"'' 


(0.5 - 1.4) X 10- 




VSGs 


Y 


(1.3 - 2.3) X 10"'' 


(1.4 - 3.8) X lO"'' 


(2.8-7.5) X 10"'' 


(1.1 - 2.8) X 10" 


-5 




a_ 


1.2 nm fixed 


1.1 - 2.5nm 


1.2 nm fixed 


1.2 nm fixed 






a+ 


3.5 - 12 nm 


2.1 - 5.0nm 


1.9 - 2.7nm 


2.0 - 3.0 nm 






a 


2.6 - 5.2 


0.2-3.7 


2.8 - 7.1 


3.2 - 36 




BGs 


Y 


(3.5 - 4.7) X 10"" 


(1.9 - 4.0) X 10"" 


(0.5 - 1.5) X 10"'' 


(0.3-1.1) X 10" 


•i 




a_ 


2.1 - 2.9 nm 


3.0 - 3.7 nm 


1.7 - 2.6 nm 


3.7- 4.0 nm 






a+ 


110 nm fixed 


110 nm fixed 


110 nm fixed 


110 nm fixed 






a 


5.8 - 35 


25-37 


4.8-15.5 


23 - 40 




VCGs 


Y 


(1.3 - 0.4) X 10"'' 


(2.0 - 0.7) X 10"^ 


(0.7-4.1) X lO-" 


(1.7-0.2) X 10" 






T 


5-7K 


7-9 K 


5-7K 


5-9 K 






P 


1 fixed 


1 fixed 


1 fixed 


1 - 2 





Table 7. Range of reliability of the dust parameters. For comparison, we give the corresponding values for NGC 1569 
(paper II). These limits take into account both the spread due to the error bars of the observations and the spread 
due to the assumed value of the radius. The indices of the power-law distribution, a, have high uncertainties since 
the size distribution is very steep. For these cases, the choice of a single grain size would effectively fit the emission. 
"Fixed" indicates that these parameters are fixed to the Galactic values due to the lack of constraints. 



of the Galaxy. In paper I we discussed the correlation be- 
tween the hardness and intensity of the ISRF and the lack 
of PAHs in the starbursting dwarf galaxies, while, PAHs 
are observed to be abundant throughout the Galaxy. We 
see an evolutionary sequence among these four galaxies. 
NGC 1140, where PAHs are evident, is the most quies- 
cent galaxy of the sample, with a radiation field that 
is less intense and relatively soft compared to the other 
galaxies (Fig.|H| paper I). In the three other dwarf galax- 
ies, NGC 1569, IIZw40 and He 2-10, PAH bands are very 
weak, if present at all (paper I). Likewise, the ISRFs are 
intense and hard. Among these three starbursting dwarf 
galaxies, NGC 1569 has a softer ISRF and the aromatic 
bands seen in the CVF spectrum are present, even if very 
weak, contrary to IIZw40, where there is no evidence for 
the PAH band emission in the CVF spectrum. 

4.3.2. The dust emission 

The best modeled dust SEDs, after iteration, are plot- 
ted in Fig. 121 and compared to that of the Galaxy in 
Fig. The characteristics of the 3 newly-presented dust 
SEDs have many similarities to those already deduced for 
NGC 1569, in paper II. 

The SEDs indicate overall hotter dust, peaking at 
^ 60 i-iTn. Due to the small grain size distributions of the 
PAH, VSG and BG components, most of the grains are 
stochastically heated and are not in thermal equilibrium 
with the radiation field in these galaxies. Table |H1 contains 
the temperature range of each component for the maxi- 
mum and minimum sizes compared to Galactic values. In 
paper II, we estimated the transition grain radius, aj, be- 
tween stochastic heating and thermal equilibrium to have 
an idea over what range of dust sizes it is important to 
consider the process of stochastic heating. This is a func- 
tion of the radiation field, the dust size and the heat ca- 
pacity of the dust. We refer to Section 4.1.3 of paper II 
for a detailed explanation of the calculation. In Fig. |H1 we 



compare the variation of the cooling rates and the photon 
absorption rates as a function of the radius, for a single 
BG. The point where the two curves intersect gives an idea 
of the grain radius, at, below which the stochastic heating 
becomes dominant. These transition radii are at ~ 5 nm in 
NGC 1569, at ~ 4nm in IIZw40 and He 2-10, Of ~ 6nm 
in NGC 1140 and at — 9nm in the Galaxy. These esti- 
mates demonstrate that the grains become stochastically 
heated at smaller radii, due to the higher intensity of the 
radiation field (Fig. O . By comparing the estimates of at 
with the typical radius that dominates the mass spectrum 
of the grains, which is a ~ 3 — 4nm (Fig. we deduce 
that most of the grains are stochastically heated in these 
galaxies. Even at wavelengths as long as 60 /im, where the 
dust emission peaks, the grains are primarily stochasti- 
cally heated. 

Another striking feature of the dust emission spec- 
tra is the weakness or absence of the mid-IR PAH 
emission bands particularly in NGC 1569 and IIZw40 
(the mid-IR spectra are presented in paper I). He2- 
10 has i iot been observed with the ISOCAM CVF. 
However. iMartin- Hernandez et all l)2005l private commu- 
nication from Marc Sauvage) have obtained a ground- 
based mid-IR spectrum in the range 8 — 13 /im which 
also reveals a lack of PAH features. Therefore, for these 
three galaxies, we only put an upper limit on the PAH 
mass, given by the continuum emission of this compo- 
nent. We do not actually fit the bands in the model. In 
the case of NGC 1140, we see relatively significant PAH 
emission, allowing us to properly constrain the PAH mass. 
For this purpose, we use a modified version of the DBP90 
model (Laurent Verstraete, private communication) where 
the PAH emission features are modeled more precisely 
(Fig. |3J). The optical constants are deduced from a spec- 
trum of a typical Hii region. For NGC 1140, we use the 
process of decomposition of the ISOCAM mid-IR spec- 
trum described in paper I: the PAH features are mod- 
eled as Lorentzian bands; the ionic lines are modeled as 
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NGC 1569 


IIZw40 


He 2-10 


NGC 1140 


Galaxy 






Tmin 


Tmax 


Tmin 


Tmax 


Tmin 


Tmax 


Tmin 


Tmax 


rmin 


Tmax 


PAH 


a_ 


2.7 K 


9200 K 


2.7 K 


9200 K 


2.7 K 


9200 K 


2.7 K 


9200 K 


2.7 K 


4400 K 




a+ 


2.7 K 


1800 K 


2.7 K 


1800 K 


2.7 K 


1800 K 


2.7 K 


1800 K 


2.7 K 


1100 K 


VSG 


a- 


2.7 K 


890 K 


2.7 K 


490 K 


2.7 K 


890 K 


2.7 K 


890 K 


2.7 K 


630 K 




a+ 


2.7 K 


160 K 


2.7 K 


350 K 


2.7 K 


380 K 


2.7 K 


340 K 


2.7 K 


78 K 


BG 


a- 


2.7 K 


230 K 


2.7 K 


140 K 


2.7 K 


190 K 


2.7 K 


111 K 


15 K 


22 K 




a+ 


28 K 


28 K 


30 K 


30 K 


31 K 


31 K 


20 K 


20 K 


17 K 


17 K 



Table 8. Temperatures of the PAH, VSG and BG components for NGC 1569 (paper II), IIZw40, He 2-10 and 
NGC 1140 compared to those of the Galaxy (DBP90). We give the minimum (Tmin) and maximum (Tmax) tem- 
peratures for each component. All the grains, except the largest BGs (for which = Tmax = Toquiubrium) , are 
stochastically heated. The values of Tmax given for the PAHs are larger than the vaporisation temperature, thus they 
are unphysical. The minimum temperature never reaches values below Tcmb, since all grain sizes are in equilibrium 
with the very low energy input coming from the CMB. 



Gaussian emission lines and the continuum is fitted as a 
modified black body providing a constraint on the VSG 
component in the dust model. Then, we add the PAH part 
of the model to fit the entire spectrum, including the aro- 
matic bands. 

For each SED of our sample, we are left with a 
submillimetrc/millimetre excess that we cannot explain 
with the standard DBP90 model. As originally found for 
NGC 1569, we are compelled to invoke the presence of an 
ubiquitous very cold grain component (VCGs) in II Zw40, 
He 2-10 and NGC 1140. Various hypotheses for this excess 
are fully explored in paper II. They include: very cold dust, 
change of grain optical properties at long wavelengths and 
grain-grain coagulation. We modeled the VCG component 
using a modified black body fSect. H!T|l . The temperatures 
of these VCG components are T = 5 - 9 K (Tabled). We 
have shown in paper II (Sect. 4.5) that the VCG compo- 
nent could correspond to very cold dust embedded in very 
dense clumps. Other means to increase the submillimetre 
emissivity which we see in these galaxies could be differ- 
ent optical properties due to temperature effects, as well 
as grain-grain coagulation. The temperature dependent 
optical properties are considered in paper II and in the 
Appendix IbI of the present paper, but so far failed to pro- 
duce better fits. Thus, in paper II , as well as the present 
paper, we chose to discuss the very cold dust hypothe- 
sis, since it has a number of quantifiable consequences 
(Sect. EH). 

4.4. The extinction curves 

From the dust properties described by the parameters in 
Tables|niand|7| we synthesize an extinction curve assuming 
a simple screen attenuation of the radiation. Construction 
of the extinction curve is based on the DBP90 model as- 
sumptions of the PAHs being the carriers of the FUV 
non-linear rise; the carbonaceao us VSGs giving rise to 
the 2175 A extinction bump fe.g. ISavaee fc Mathis 197^) 
and the BGs giving rise to the NIR and visible rise 
of the extinction curve. The individual contributions to 
the extinction curves and the total synthesized extinc- 



tion curves are displayed in Fig. El and are compared 
to the Milky Way (DBP90 ), the Large Mag ellanic Cloud 
l|Koornneef fc Coda.l981t.N'aiidv et al.lll98ll) and one line 
of sig ht toward the Small Magellanic Cloud l)Prevot et alJ 
Il984l) in Fig.Uni 

These synthesized global extinction curves are char- 
acterised by a quasi-linearity with 1/A which is also the 
case for the extinction curves seen toward the Magellanic 
Clouds fFig. irnj) . These slopes indicate a higher absorp- 
tion in the UV as a consequence of the presence of the 
small sizes of the BGs, as can be seen on the decom- 
positions in Fig. El and also due to the lack of PAHs 
which are responsible for the far-UV rise. We also note 
that the 2175 A extinction bump is lower in NGC 1569, 
He 2-10 and NGC 1140 than in the Milky Way, as is the 
case for the extinction curves observed toward a line of 
sight in the low-metallicity Magellanic Clouds. This is 
a consequence of the higher silicate-over-graphite ratio 
(Mbg/A^vsg); in these galaxies than in the Milky Way 
(this ratio is equal to 14 for the Milky Way, and to 24, 18 
and 28 for NGC 1569, He 2-10 and NGC 1140, respectively, 
from Table El). Indeed, in the DBP90 model, the carriers 
of the bump are the VSGs. However, the ratio Mbg/Mvsg 
is smaller in IIZw40 (A/bgA^vsg — ll)j which explains 
why the model determines a larger bump in IIZw40. The 
silicate-graphite ratio is not constrained very accurately 
as can be seen from the error bars on our parameters 
in Table [7\ Moreover, the carriers of the bump could be 
other specie s of g rains like the PAHs, as proposed by 
Joblin ct al. {1992). In recent dust mode ls l|Draine fc Li 
2001: Li fc Draine.,200lllZubko et alJl2004j) . the carriers of 
the bump are both PAHs and solid carbonaceous grains. 
Thus, the lack of PAHs could naturally explain the weak- 
ness of the bump in these galaxies. 

We think, a s already suggested e.g. by 
iMas-Hesse fc KunthI |T999) and IWeingartner fc Drain j 
l)200ll) . that the shape of the extinction curve, thus the 
very different dust properties, essentially reflects the 
starburst activity in these galaxies, rather than a conse- 
quence of their low metallicity. Indeed, the fragmentation 
and the erosion of the grains by the numerous shock 
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Fig. 3. IIZw40, He 2-10 and NGC 1140 observations and 
modeled SED. The data (TableQ are indicated by crosses: 
vertical bars are the errors on the flux values and the hor- 
izontal bars indicate the widths of the broadbands. The 
lines are the dust model and the different dust compo- 
nents. Diamonds are the model integrated over the ob- 
servational broadbands and color-corrected. Thus, depar- 
tures from the model lie where the diamonds deviate from 
the crosses. The power is expressed in WH~^ which is i^L^ 
divided by the number of H atoms (L^ is the monochro- 
matic luminosity). 



Fig. 4. Size distribution deduced for IIZw40, He 2-10 and 
NGC 1140. The solid thick lines are the best fit solution 
and the non-solid lines are the solutions for the two ex- 
treme values of the radius. The grey thin line is the size 
distribution for our Galaxy (DBP90). 



waves following the starburst, as suggested here, could be 
the main process determining the extinction properties, 
rather than a deficit in formation of large grains by 
evolved stars. This argument is also supported by the 
fact that observations toward some lines of sight in the 
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SMC show Galactic type extinction ijGordon &: Clavtonl 
Il998t lOordon et al J 12003) and observations toward some 
lines of sight inside the Gala xy show SMC like extinction 
curves l|Valencic et al.lEo03|) . 

4.5. The dust masses 

Table El contains all of the most interesting physical quan- 
tities deduced from our modelling. The dust masses are 
directly deduced from Tables [HI and [71 Table highlights 
the fact that most of the mass is contained in the VCG 
component, if we believe that the submillimetre excess is 
produced by very cold dust. There is a large uncertainty 
on the mass of this component since its temperature is 
uncertain. However, our results show that between 40 and 
80 % of the dust mass is in the form of VCGs. The rest of 
the dust mass is dominated by the BGs. The VSGs, while 
the most numerous particles, contain between 1 and 3 % of 
the total dust mass. To compute the gas-to-dust mass ra- 
tio we assume that Mgas = M (H i) + M(H2) + M(Re) and 
that M(He) = 0.25 x Mgas- The M(Hi) are taken from 
tableland Af (H2) are given in the literature: iMeier et alJ 
[2001,,) report M(H2) < 0.4 x 10^ Mq for nZw40 and 
M(H2) = 1.4 X 10* Mq for He 2-10. For NGC 1140, we as- 
sume that A/(H2) is negligible compared to A/ (Hi), since 
we did not find any global estimate in the literature and 
this mass is likely to be small. We use the following defi- 
nition of the dust-to-metal mass ratio: V = l/QZ, which 
is the fraction of metals locked up in dust. The Galactic 
values are typically ~ 150 and V ~ 1/3. The range of 
values we give for Q and V are rather large (Table 0. The 
constraints that we put on this parameter allows lower 
values of the du st-to-metal mass ratio th an established 
by other studies ijLisenfeld fc Ferraralll998j) . The dust-to- 
metal ratio for the dwarf galaxies may be lower than that 
of the Milky Way, except for the upper limit of II Zw 40. It 
would suggest that low-metallicity systems could be less 
efficient in forming dust than normal metallicity galaxies. 

He 2-10 is a special case, since the VCG mass we com- 
puted is affected by the fact that the upper limit of the 
free- free emission we subtracted from the millimetre fluxes 
in order to be conservative, is significant, and also, since 
its metallicity is relatively uncertain. This is the reason 
why the value of V in He 2-10 deviates from the others. 

4.6. The nature of the very cold grain component 
(VCGs) 

In paper II, we discussed the likeliness of the VCG hypoth- 
esis by estimating order of magnitudes of the temperature 
and the filling factor of dust hidden in very dense clumps. 
Here, we do the same calculation for II Zw40, He 2-10 and 
NGC 1140. 

Following the equation (11), in paper II, we compute 
a rough estimate of the temperatures of the VCGs, TycG- 
This is the temperature that would be reached by grains 
totally shielded from the stellar radiation, heated only by 



the emission from the surrounding dust. It is estimated 
assuming that the inferred mass of VCGs will be heated 
only by the radiation arising from the dust peaking in the 
FIR in our synthesised SEDs. The ranges of values that 
we calculate are: Tyco — 7 — 9K for IIZw40, TycG — 
8 - lOK for He 2-10 and TycG ~ 6 - 9K for NGC 1140. 
All these values are consistent with the temperatures given 
in Table except in the case of He 2-10 where the high 
level of the upper limit of the free-free emission gives more 
uncertainty on the VCG parameters. We believe that these 
figures support the hypothesis of very cold dust in dense 
clumps to explain the submm excess of the SEDs. 

We are also able to give an idea of the size of the 
clumps and of their filling factor. As we did in paper II, 
for NGC 1569, we assume that the BGs are shielding the 
VCGs. If we assume typical contrast densities oiS = 10"* — 
10^, we find clump sizes of 0.9 - 18 pc for IIZw40, 1.2 - 
38 pc for He 2-10 and 2.4 - 43 pc for NGC 1140. The corre- 
sponding fiUing factors are: 2 x 10^** — 2 x 10^"^ for II Zw 40, 
3 X 10"'^ - 3 X IQ-^ for He 2-10 and 4 x 10"^ - 8 x 10""^ for 
NGC 1140. Thus, our conclusion is that the ISM of these 
galaxies could be very clumpy with low filling factors. In 
contrast, the filling factor of the Gala ctic molecular phas e 
is about 1 % l|Tielen£lll995|) . Recentlv. lA"ndre et all l)2004j) . 
found very small (~ 100 pc) dense CO clumps in the low- 
metallicity ISM of the Magellanic clouds. These clumps 
could cont ain the very cold du st probed by our submm- 
mm SEDs. iRubio et all l|2004|) found that the CO mass 
deduced from the submillimeter continuum emission of 
molecular cloud in the SMC, is ~ 10 times higher than 
the virial mass reported from CO observations. Their in- 
terpretation is that the CO emission comes from dense 
clumps. 

5. Summary and conclusion 

We have presented new SCUBA images at 450 /im and 
850 Aim and MAMBO ON-OFF observations of IIZw40, 
He 2-10 and NGC 1140. With additional data from the ht- 
erature, we have constructed the observed SEDs for these 
galaxies. Using the same modelling approach as in pa- 
per II, we have combined stellar evolution, photoionisa- 
tion and dust modelling to compute self-consistent dust 
SEDs. The results found in paper II, for NGC 1569, are 
supported here, for IIZw40, He 2-10 and NGC 1140. 

1. We find very low abundances of PAHs, and smaller 
overall sizes (~ 3 — 4nm) of grains emitting in the 
mid-IR and far-IR, in contrast to grains observed in 
the Galaxy. The small sizes of the grains, on average, 
are supported by the erosion by shock waves produced 
by the numerous supernovae that occured in these star- 
burst galaxies. 

2. Due to the small sizes of the grains, the stochastic 
heating of dust is predominant in these galaxies, even 
the grains emitting at far-IR wavelengths. The bulk 
of the emission comes from grains which are not in 
thermal equilibrium with the radiation field. 
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NGC 1569 


II Zw40 


He 2-10 


NGC 1140 


Ldvst (L©) 
LT' (Lo) 


(5.8 ± 0.1) X 10" 
(1.3 ±0.1) X 10^ 
0.24 


(2.0 ± 0.2) X 10" 
(2.5 ±0.1) X 10^ 
0.15 


(5.4 ±0.2) X 10" 
(5.3 ± 0.2) X 10^ 
0.26 


(3.3 ± 0.1) X 10" 
(6.0 ±0.2) X 10^ 
0.11 


MpAH (M0) 
MvSG (M0) 
Mbg (Mq) 

MvCG (M0) 

Total Md (Mq) 


^ 190 

3.4lJ° X 10^ 
8.4l?--^ X 10* 
(0.7-2.5) X l(f 
(1.6 - 3.4) X lO'' 


< 220 

1.0t°;I X 10* 
l.lt^ol X 10^ 
(3.1-8.8) X lO'' 
(0.4 - 1.1) X lO'' 


< 930 

1.5l°;«7 X 10* 
2.%t\% X 10^ 
(0.2 - 1.3) X lO*' 
(0.4 - 1.8) X 10® 


6.2+^-^ X 10^ 
1.7tii X 10« 

A.etli X 10® 

(0.2 - 1.3) X lO'' 
(0.5 - 2.0) X lO'^ 


d 

V 


740 - 1600 

1/6 - 1/3 


530 - 1460 

1/5 - 1/2 


330 - 1500 

1/30 - 1/6 


500 - 2000 
1/13 - 1/3 



Table 9. Derived quantities from our dust modelling, for NGC 1569 (paper II), IIZw40, He 2-10 and NGC 1140. Ldust 
is the luminosity reemitted by the dust from near-IR to mm. LJ'^* is the escaping stellar luminosity. Ty^ is the effective 

optical depth deduced from the energy balance (sec paper II). AfpAH, ^'-^VSG, Mbg and AfycG arc the masses of each 
dust component and Mj, is the total dust mass. G is the gas-to-dust mass ratio and V is the dust-to-metal mass ratio. 
The Galactic values are typically Q ~ 150 and 2? ~ 1/3. 



3. In each of our four SEDs presented here and in pa- 
per II, we find, with no exception, a submillimetre 
emission excess. Wc propose that this excess is the 
emission of very cold dust (5 K < T < 9 K) hidden in 
dense clumps. In each of these SEDs, except one, we 
find that P = 1 is the most likely value for the emis- 
sivity index of this component, even if this parameter 
is difficult to accurately constrain with the scarcity of 
data in the submillimetre-millimetre part of the spec- 
tra. Very cold grains represent between 40 and 80 % 
of the total dust mass. We show that this very cold 
dust hypothesis is consistent, although we can not ex- 
clude other explanations for this millimetre excess, like 
non-standard optical properties. 

4. The submillimetre-millimetre emission is not concen- 
trated toward the outer regions of the galaxies but cor- 
related with the mid-IR emission distributed around 
the star forming regions. Our results are consistent 
with a clumpy medium composed of ubiquitous clumps 
of sizes between a few pc to a few tens of pc with typ- 
ical filling factor of lO"" - lO'^. 

5. The extinction curves that we synthesize for the 4 
dwarf galaxies are different from that of the Galaxy 
and similar in shape to that of the LMC. Due to 
the small sizes of the grains and the low abundance 
of PAHs (or lack of PAHs), the slope of the extinc- 
tion curves are constant from near-IR to far-UV wave- 
lengths with the greatest extinction in the UV range. 
The extinction bump at 2175 A, which is assumed to 
be carried by the very small grains in the dust model 
we use, is found to be smaller than the Galactic one 
in our sample except in one case where we have a high 
mid-IR attenuation. 

6. We constrain the gas-to-dust mass ratios which lie be- 
tween 300 and 2000 and the dust-to-metal mass ratios 
which lie between 1/30 and 1/2 in our sample. 

This paper closes a serie of three publications related 
to the MIR to mm properties of low-metallicity environ- 
ments. This study investigates the dust properties of 3 new 



low-metallicity galaxies, IIZw40, He 2-10 and NGC 1140 
and along with NGC 1569 (paper II), presents the first de- 
tailed SED studies focused on dwarf galaxies. The study 
of dust properties using IR emission provides results that 
are overlooked by optical studies based on the extinction 
phenomenon. Dust properties are profoundly different in 
these galaxies compared to our Galaxy and other nor- 
mal galaxies. The dust size distribution is different, con- 
sequently the SEDs and the extinction laws are not the 
same. Moreover the dust mass and temperature estimates 
arc affected by both the presence of very cold dust and 
by the non-thermal equilibrium heating of the bulk of the 
emitting dust. Assuming Galactic dust properties when 
looking at these galaxies is not correct and could lead to 
large errors, particularly in dust masses. 

We also warn the reader that it is very important to ob- 
tain data in the MIR regime as well as in the submillime- 
tre/millimetre regime to properly constrain the dust prop- 
erties. The submillimetre/millimetre regime is still largely 
unknown and poorly sampled. Observations of galaxies 
at these wavelengths could contribute significantly to the 
comprehension of the variety of galactic dust environ- 
ments. Without such precise dust SED modeling incorpo- 
rating the submillimetre/millimetre wavelength range, in- 
accurate results could lead to poor assumptions in galaxy 
number counts and galaxy evolution models. The future 
ground-based, airborne and space observatories exploring 
the FIR to millimetre wavelength regime, such as SOFIA, 
ASTRO-F, Herschel, Planck and ALMA hold great expec- 
tations in the investigation of extragalactic dust properties 
and will bring us many steps closer ot the understanding 
of the physics of the ISM of galaxies. 
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Appendix A: Degeneracies in dust modelling 

In this appendix, wc demonstrate a well-known degener- 
acy between the energy source, and the grain properties, 
which must be considered when modelling a dust SED. 

The dust SED of a given region depends (i) on the 
optical properties of the different types of grains, (ii) on 
the size distribution of these grains and (Hi) on the shape 
and the intensity of the heating radiation field. 

Optical properties of the grains are based on labora- 
tory experiments and astrophysical measurements, it is 
not relevant to vary them in order to fit an observed SED. 
However, the radiation field depends on the environment. 
The grains will experience a different radiation field in a 
starburst region than in the diffuse ISM. The top-left plot 
of Fig. IX!n show s the fit of a fake ob served SED (the grey 
crosses) with the lDesert et all ()l99Cl|) model by varying the 
ISRF but keeping the Galactic grain size distribution. The 
radiation field needed to fit this SED is the grey line on 
the bottom-left plot on Fig. lA.ll compared to the Galactic 
one, in black. 

The dust size distribution is also a set of parameters 
varying significantly in different environments. Indeed, it 
is subject to fragmentation and erosion by shock waves, 
evaporation near strong radiation sources, or coagulation 
and accretion of material in dense media. The top-right 
plot of Fig. lA.ll shows the same fake o bserved SED than on 
the top- left plot, but fitted with the iDesert et all l)l99Cl|) 
model by varying the size distribution and keeping the 
Galactic radiation field. We get a good fit, in this case, 
too. The bottom-right plot of Fig. IXTl shows the required 
size distribution to fit the SED, in grey, compared to the 
Galactic one, in black. 

These two examples show that, without caution, we 
can interpret the same SED with two different physical 
solutions. We emphasize the need to constrain both the 
ISRF and the size distribution when modelling a dust 
SED. Our modelling approach was aimed to achieve that 
goal. 

Appendix B: Peculiar optical properties to explain 
the submillimeter excess 

lAgladze et al.l lll996() measured the variation of the emis- 
sivity index, f3, of various silicate and carbonaceous ma- 
terials with the temperature (Fig. IB.1|) . We have imple- 
mented their observations in order to explain the submil- 
limeter excess of our SEDs. 



We modified the optical properties of the BGs, by con- 
structing a temperature dependent absorption efficiency 
Qabs(A,a,T) (Fig. |R2l). In the standard DBP90 model, 
the BGs are silicate grains and have an index (3 = 2, in 
the submillimeter regime. The amorphou s MrSIOs and 
Mg2Si04 are the two species studied by lAgladze et alJ 
which are consistent with the BGs. This modifi- 
cation can have a significant effect, if the contribution by 
cold BGs (Tbg < 25 K) is important. 

For a given grain size, cold BGs can exist due to tem- 
perature fluctuations. They cool down between two pho- 
ton absorptions. In our case, the monochromatic luminos- 
ity emitted by a BG of radius a is: 

/■°° dP 
Cx{X,a)= / ATTBx{X,T)Tra^Q,^,{\,a,T)—dT, 

(B.l) 

where B\{X,T) is the Planck function, and dP/ dT the 
temperature distribution. Otherwise, BGs can be cold if 
they are large enough. The total luminosity emitted by 
the BGs, with the number size distribution /(a), is: 

/>oo 

Lx{\) = / Cx{X,a)f{a)da. (B.2) 
Jo 

Among our four dwarf galaxies, the only detectable 
effect of these peculiar optical properties can be seen in 
NGC 1569 (Fig.|R3l). For the three other SEDs, the two 
models give identical results. 

Finally, we conclude that this effect is not able to re- 
produce the submillimeter excess. This does not mean that 
we can exclude other peculiar, as yet unexplored, optical 
properties to explain this excess. 
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Fig. 6. Comparison of the synthesized ISRFs of NGC 1569 
(paper II), IIZw40, He2-10 and NGC 1140, with the 
Galactic one (DBP90). These ISRFs are the effective ra- 
diation fields seen by the dust in our model. 
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Fig. 7. Comparison of the modeled dust SEDs of 
NGC 1569 (paper II), IIZw40, He2-10 and NGC 1140 
(this paper), and that of the Galaxy (DBP90). 
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Fig. 5. Synthesized ISRFs for IIZw40, He 2-10 and 
NGC 1140 computed with PEGASE and CLOUDY. The 
points with error bars are the observational data from 
Table O the solid black line is the global non-extincted 
ISRF, the dashed line is the global extincted ISRF and 
the grey line is the young single-burst component. The 
extinction curve used is the output from the dust model 
DBP90. 
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Fig. 8. Cooling rates and photon absorption rates as a 
function of the grain radius a for NGC 1569, II Zw 40, He 2- 
10 and NGC 1140, compared to the Galaxy. 
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Fig. 9. Synthesized extinction curves for llZw40, He 2-10 
and NGC 1140. The contributions to the extinction are 
shown individually for the PAH, VSG and BG compo- 
nents. The solid lines are the total extinction curves mod- 
eled with DBP90. The opacity is expressed as the cross 
section per H atom. 
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Fig. 10. Comparison of the synthesized extinction curves Fig. B.2. Modification of the opti cal properties of the 
of NGC 1569 (paper II), IIZw40, He2-10 and NGC 1140 BGs, following the measurements bv lAgladze et alJ jli 
with the Galactic extinction curve (DBP90). the ob- shown in Fig. IB. II 
serv ed LMC curv e (aver age of iKoornneef fc Cod i lll98lh 

and iNandv et all ' 198l|) 'l and the observed SMC curve NGC 1569 

i Prevot et al.lll984 1. These curves have been shifted for 
clarity. 
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Fig. B.3. Fit of the long wavelength SED of NGC 1569 
(paper II), with the standard model (in black), and with 
the modified optical properties (in grey), using the tem- 
perature dependent Qabs- 



Fig. B.l. Variation of the emissivity index, /3, as a func- 
tion of the temperature, T, for Mg2Si04, as measured by 
fAgladze et al. (,1996.) . 



